Detection of qubit-oscillator entanglement in nanoelectromechanical systems 
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Experiments over the past years have demonstrated that it is possible to bring nanomechanical 
resonators and superconducting qubits close to the quantum regime and to measure their properties 
with an accuracy close to the Heisenberg uncertainty limit. Therefore, it is just a question of time 
before we will routinely see true quantum effects in nanomechanical systems. One of the hallmarks 
of quantum mechanics is the existence of entangled states. We propose a realistic scenario making 
it possible to detect entanglement of a mechanical resonator and a qubit in a nanoelectromechanical 
setup. The detection scheme involves only standard current and noise measurements of an atomic 
point contact coupled to an oscillator and a qubit. This setup could allow for the first observation 
of entanglement between a continuous and a discrete quantum system in the solid state. 

PACS numbers: 85.85.+j, 03.67.Mn, 72.70.+m 
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In recent years, nanoelectromechanical systems 
(NEMS) have become a strong focus of research in the- 
oretical and experimental physicsji One of the practical 
reasons for this development is the prospective use of 
NEMS to design devices which allow the measurement of 
position, force and masa^^— with unprecedented accura- 
cies. From a more fundamental point of view, NEMS op- 
erate at the boundary between the classical and the quan- 
tum world and recent works have suggested that NEMS 
will soon allow the observation of quantum mechanical 
states in mesoscopic mechanical systems.— 

The observation of quantum states of matter in such 
systems generally requires ultralow temperature and low 
dissipation. Bringing a nanomechanical oscillator of fre- 
quency f2 near its ground state means reaching a temper- 
ature T ^ hO. Iks- Various schemes to cool an oscillator 
to its ground state have been proposed;^^— and experi- 
ments on nanomechanical systems are now approaching 
this limit. ^'^^ Moreover, high quality factors have been 
achieved which lead to relaxation and decoherence times 
long enough for the measurement of quantum states.— 

One of the most rewarding endeavors involves the cre- 
ation and detection of nonclassical correlations (entan- 
glement) between the nanomechanical oscillator and an- 
other quantum system. The easiest option would be to 
entangle the oscillator with a mesoscopic system whose 
properties are well understood and in which quantum 
effects can routinely be observed: a superconducting 
qubit. It has been demonstrated that these devices 
have decoherence times which can exceed oscillation pe- 
riods of nanomechanical resonators by several orders of 
magnitude!^ Various theoretical proposals have been 
made on how entanglement between an oscillator and 
a qubit can be createcU^^ and such systems have been 
successfully coupled in experiments.— 

In this article, we propose a system which allows the 
detection of entanglement between an oscillator and a 
qubit using an electronic measurement in an atomic point 



contact (APC). The electronic system is based on a tun- 
neling contact, a readout device which is known to be 
quantum-limited.— We find that the measurement of the 
current and the symmetrized current noise in this system 
allows the evaluation of a criterion for entanglement^ 
based on the density matrix of the oscillator-qubit sys- 
tem. This allows for the detection of entanglement in 
arbitrary pure or mixed states. All elements of the pro- 
posed setup have been realized separately in different ex- 
periments. Moreover, it has been shown that the current 
and the noise of an APC can be measured with a high 
accuracy. Therefore, it should be possible to combine 
both elements into one functional device as schemati- 
cally shown in Fig. [T] and to measure its current and 
noise properties. 

The system we investigate consists of a nanomechan- 
ical oscillator, a qubit and a biased APC. Both the os- 
cillator and the qubit are coupled to the APC and thus 
modulate its transmission coefficient. The APC consists 
of two electron reservoirs ( "left" and "right" ) at chemical 
potentials [Il.r which are subject to a voltage difference 
V = [i]^ — [Lji. The Hamiltonian of the APC reads (using 
units e = fi, = 1) 

a=R,L k 

Ht = lY^pliz = 0)4,r{z - 0) + h.c. , (1) 

where i^ti^) = X]fc ^^^'^^V'l fc — L,R) creates 

an electron at position z in the left/right reservoir, re- 
spectively. As a simplification, we assume a constant 
density of states po = I/Ittwf) where uj? is the Fermi 
velocity. In the tunneling Hamiltonian Ht, the count- 
ing operator Y (Y^) decreases (increases) the transferred 
charge by one, Y \nfi + \) = |ni{). The corresponding 
number operator is defined by fin jriij) = \niij. The 
tunneling amplitude 7 will be specified shortly. 

The oscillator and the qubit are described by the 
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Hamiltonians 
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Hq 



(2) 



where m and $7 denote the effective mass and the fre- 
quency of the osciUator, respectively. In the qubit Haniil- 
tonian, cr^^y^z denote the Pauh matrices. For A = 0, the 
energy difference between the two qubit states is given 
by 2e. A finite A enables tunneling between the states. 

The state of the qubit-oscillator system modulates the 
tunneling amplitude 7 of the APC. If the oscillator acts 
as one of the electron reservoirs of the APC^^ as shown 
in Fig. [TJ the tunneling gap depends on the oscillator 
displacement x. For small x one obtains 7 oc 70 + 71a;. 
The same dependence can also be realized for capacitive 
coupling f2£ The qubit can be realized as a Cooper pair 
box in which case a depletion of the electron reservoirs 
of the APC depending on the state of the qubit leads 
to an additional term ^2<^z in the tunneling amplitude. 
Irrespective of the concrete realization, to lowest order 
the combined effect of the oscillator and the qubit leads 
to 



7 = 70 + 712; + 720-Z 



(3) 



In general, the amplitudes jj — |7j|e"^^ (j — 0,1,2) can 
be complex. Since the global phase is irrelevant, we set 
i5o = 0. Finite phases 61^2 can be realized experimentally 
by closing the electric circuit using an additional tunnel 
junction as shown in Fig. [I}S^ Threading the loop with a 
magnetic flux causes Aharonov-Bohm phases which can 
be absorbed in the tunneling amplitudes and generally 
lead to finite phases Si and S2 ■ This is discussed in more 
detail in Appendix |^ The benefits of a controllable Si 
have been investigated for a system consisting of an APC 
and an oscillator i^"^ while for Si — 0, the current noise 



only depends on the oscillator position 



a finite Si 



leads to terms proportional to (p'^) and thus contains 
information about the oscillator momentum. Similarly, 
the presence of tunable phases Si,2 increases the number 
of measurable oscillator and qubit properties. 

The current operator is defined by / = 
main objective will be the calculation of the average cur- 
rent (/) and the symmetrized noise spectral density. 
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dTe-(({/(r),/(0)})) 



(4) 



where ((/(r)/(0))> = (/(r)/(0)) - (/(r)) (/(O)) and {•,•} 
denotes the anticommutator. We shall find that the av- 
erage current as well as the noise depend on expectation 
values of products of qubit and oscillator operators. From 
these expectation values, the bipartite expectation value 
matrix (EVM) can be constructed. It is a complex 6x6- 
matrix defined by^*' 
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(5) 




FIG. 1: (Color online) Possible experimental setup consist- 
ing of a qubit and an oscillator coupled to an atomic point 
contact (APC). Electrons tunnel at the APC (ta) and a fixed 
tunnel junction (tb), which are both biased with a voltage 
V. The area enclosed by the junctions (red dashed line) is 
threaded with a magnetic flux to create an Aharonov-Bohm 
phase. The qubit is realized as a Cooper pair box (CPB, yel- 
low). Its state can be tuned using the gate voltage Vg and it 
couples capacitively to both junctions. The oscillation of the 
nanomechanical resonator (NR, green) modulates the tunnel- 
ing amplitude ta- As discussed in more detail in Appendix 1X1 
this setup can be used to realize the tunneling amplitude ((3|. 



where 
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and we used Sxp = ^{xp + px). All expectation val- 
ues are taken with respect to the qubit-oscillator state 
described by the density matrix Pq^osc, e.g. (|t) (tl a;) = 
Tr[pq_osc It) (tl x]- It bas been shown^° that for any sep- 
arable state Pq,osc, 








-1 > 







(7) 



i.e. both matrices must be positive semidefinite. Here, 
Pq = Trosc Pq,osc dcuotcs the reduced density matrix of 
the qubit. Once the complete EVM is known, a violation 
of Eq. (O proves that Pq^osc is an entangled state. In the 
following, we show that the current and its noise contain 
enough information to construct the EVM and thus de- 
duce entanglement of the oscillator and the qubit. Hence, 
this measurement provides a separability criterion for an 
oscillator and a qubit, comparable to a Bell inequality 
measurement which provides a separability criterion for 
two qubits. 

Since the oscillator-qubit state is not necessarily sta- 
tionary, both the average current and the noise spec- 
trum will in general be time dependent. Therefore, we 
relate the current and noise at time t to the expecta- 
tion values of the qubit-oscillator system taken at the 
same time. We calculate the current using perturba- 
tion theory in the tunneling Hamiltonian Ht and using 
Hq = Hei + -fffl -l- Hq as the unperturbed Hamiltonian. 
The method is discussed in more detail in Appendix |B] 

Calculating the time-dependence of the qubit and os- 
cillator operators using the Hamiltonian Hq, the Kubo 



formula straightforwardly yields the average current, 



(/) = 2VT^ + -^r2 + -^Ti if) 



- y/JhTi [2V cos((5i ){x}+n sin(Ji ) (p) ] 

- 2 v/ll^T^ [V cos((52) ((7,) + A sin(52) {(Ty)] 
'"^ Fcos((52i) (icr^) - §sin((52i) (pcr^ 



- A sin((52i 



(8) 



where we used Six = S2 ~ Si. The dimensionless con- 
ductance of the bare tunnel junction is given by Tq = 
^PoItoP- Moreover, using the relaxation rates for the 
oscillator and the qubit^^i^i 



ri-27rp2|^i|2/m, 
T2 ^ inplVl-fif , 



(9) 



we defined dimensionless conductances for the oscilla- 
tor and the qubit by Ti = Ti/fl and T2 = T2/V. 
The oscillator displacement is measured in units of its 
zero-point motion xq = l/VSmfi, i.e. x 
p = 2xop. Finally, /q.osc = 



x/xo and 
iriSgn(T/)min(|T/|,17) is 



the current due to the zero-point fluctuations of the os- 
cillator. A similar term arises for the qubit, /o.qd^l > 
2a) = T2A {a,) sgn{V) whi le /o,qub it(|^| < ' 2a) = 
T2V (a^) A/(2a) where a = + A^. 

For 72 = 0, the result coincides with the known re- 
sult for a system containing only the oscillator coupled 
to an APCr^ while for 71 = (52 = 0, the known result 
for an APC coupled to a qubit emerges <^ In the gen- 
eral case, a current measurement enables us to deduce 
the qubit expectation values needed for the construction 
of the reduced density matrix pq. However, the correla- 
tion functions contained in the result for the current are 
insufficient to construct the EVM. 

The calculation of the noise spectral density can most 
easily be accomplished by using the Redfield equation 
in connection with the Born-Markov approximation. 
For this purpose, we split the complete system into a 
fermionic "bath", Hb = H^i and the actual "system", 
Hs — Hq + Hr, coupled by Ht- The bath timescale can 
be estimated to be of order 1/V since eV/{2-Kh) corre- 
sponds to the attempt frequency at which electrons ar- 
rive at the tunnel junction. In the limit y 3> ri. A, e, we 
can assume the bath timescale to be much faster than the 
system timescales. To lowest nonvanishing order in the 
tunneling, the equation of motion for the reduced system 
density matrix ps(t) ~ Tr^ p{t) then reads 

Ps{t) = -i[Hs,ps{t)] 

dTTYB[HT,[HT{-T),ps{t)^PB]] ■ (10) 

Jo 

This equation can be used to calculate expectation values 
of system operators. In particular, the average current 
can be calculated using (J) = ^ (nn) = Tislniips]- The 
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FIG. 2: (Color online) Schematic density plot of the prefactors 
Sx{i^) {X = xaxyXay, . . .) of the frequency-dependent noise 
^i'^) ~ X^x Sx{ijj) {X) as a function of 5i and cj. 



current calculated in this way reproduces Eq. which 
was derived without using the Born-Markov approxima- 
tion, as soon as V > 2\/€^ -I- A^. 

In order to use Eq. (|10p to calculate the symmetrized 
noise, we make use of MacDonald's formulai^^ 



S(uj) = , 



dt smiu;t)-{{nlit))) . (11) 



The calculation of the cumulant derivative -^{{hj^{t))) us- 
ing Eq. (jlOp leads to expressions containing cumulants 
{{hjiX)), where X is in general a product of system op- 
erators. The time evolution of these cumulants will lead 
to ever higher-order cumulants, such that the resulting 
set of differential equations is not closed. In order to 
make the problem tractable, several approximations are 
needed. 

Since the qubit and the oscillator are only weakly cou- 
pled to the APC, we can assume 71X0,72 ^ 7o- In this 
case, we can ignore higher-order cumulants of the form 
{{fiRx'^)), {{fiRXCTj)), etc. Setting these higher-order cumu- 
lants to zero is often referred to as the Gaussian approx- 
imation. 

Second, the Redfield equation leads to couplings be- 
tween qubit and oscillator; e.g., -^{{fijix)) contains a term 
depending on {{hRaz))- However, since such crossterms 
are subleading in terms of tunneling amplitudes, we can 
neglect them and thereby decouple the set of differential 
equations into qubit and oscillator parts. 

Finally, we set e = and S2 — 0. This is not a crucial 
approximation but it simplifies the differential equations 
considerably. 

It turns out that this set of approximations makes a 
solution of the Redfield equation possible and the noise 
can be calculated using Eq. ([TT|) . For details, we refer to 
Appendix [Cj The resulting expression S{uj) contains all 
entries of the EVM with the exception of (jp'cjx) , (jP'cjy) 
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and (SxpCTx) ■ The frequency-dependence is characterized 
by the following functions, 



pie, consider the qubit-oscillator state 



ai{uj) = 
02(0;) = 
/3i(w) = 
/32(c^) = 



Tiw^ sin((5i) - 2y(cj^ - cos(5i) 
2VTirj cos((5i) + 17(^2 _ n^) sin(5i) 



F(4A2 



r2w2 ^ (^2 ^ 4^2)2 

r2A2 

Tluj^ + (a;2 - 4A2)2 

r2 + c.2 ' 



(12) 



which contain Lorentz and Fano shaped resonances at the 
characteristic frequencies of the system, uj — 0, il, 2A. 
The complete expression for the noise reads S^cd) — 
J2x i'^) i-^) where X denotes all combinations of 
qubit and oscillator operators contained in the EVM ([5]) . 
As mentioned above, it turns out that all except three 
of the prefactors S'x(w) are nonvanishing and are dis- 
tinguishable combinations of the functions q;i^2(w) and 
/3i,2,3('j-'). Results for Sxioj) can be found in Appendix 
[Ul A plot of the relevant cross-correlations' prefactors is 
shown in Fig. [5] Since the shapes of these functions are 
rather distinct, the expectation values constituting the 
EVM can be recovered from the total measurable noise 

Beyond the regime V ^ ^1, A, the Born-Markov ap- 
proximation fails and we have to resort to conventional 
perturbation theory in order to calculate S{uj). For the 
case of a simple tunneling junction, it is well known that 
S{uj) has kinks at |w| = In the presence of the 

oscillator and the qubit, we find that S{uj) shows similar 
features at |a;| = |y ±2A| and |a;| = |y ±rj|. The ampli- 
tudes of these kinks yield additional information about 
the expectation values (xaj) and {p(7j) {j = x,y,z). De- 
tails are given in Appendix [Cl 

Even under ideal circumstances, not all coefficients of 
the EVM can be determined by a current and noise mea- 
surement. However, it was shown that even incomplete 
knowledge of the EVM allows for the detection of entan- 
glement for experimentally relevant states^S If the EVM 
contains a few unknown parameters a^, we can only de- 
tect entanglement for states where Eq. ([7]) is violated for 
arbitrary a^. Because of the special structure of the in- 
equalities d?]) , this numerical problem reduces to a convex 
optimization problem which can be efficiently solved us- 
ing semidefinite programming as shown in [2C| . A route 
to obtain all matrix elements of the EVM would be to 
rotate the qubit state by applying 7r/2-pulses using a tun- 
able gate voltage i^ii^ii^ 

The detection of entanglement based on the EVM is 
highly versatile since it works even for arbitrary mixed 
states. The problem is of course greatly simplified if prop- 
erties of the state to be measured are known, e.g. due 
to a tailored preparation of the system. As an exam- 



— po 11, n -I- 1), where \n) denotes a Fock state of the 
oscillator and < Po < 1 ■ If the system is assumed to be 
in a state of this form, detection of entanglement reduces 
to a measurement oi pq. For this state, it turns out that 
S{n, Si) — S{fl, — (5i) ^ ^po(l - Po) sin(Ji). A measure- 
ment of the noise at the oscillator resonance frequency for 
finite Si is then sufficient to detect entanglement. Simi- 
lar relations can be derived for other states. The less is 
known about the state, the more information has to be 
gained from S{uj, Si). We found that even for completely 
arbitrary states, entanglement can be detected by this 
electronic measurement. Therefore, it is our expectation 
that this setup is ideally suited to detect entanglement 
between a continuous and a discrete quantum system. 
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Appendix A: Tunneling phases 

Phases of the tunneling amplitudes can be created by 
adding a second tunneling junction to the system and 
threading the resulting loop with a magnetic flux $. If 
the setup of Fig. [1] is chosen, one of the tunneling ampli- 
tudes will be cc-dependent while the other one will not. A 
schematic drawing is shown in Fig. |3^, and the tunneling 
amplitudes of the two junctions are given by 



ta — taO + talX + ta2<7zi 



(Al) 



where x = x/xq. Due to the Aharonov-Bohm effect, 
the magnetic flux leads to a phase shift in the electronic 
wave function which can be absorbed in the tunneling 
amplitudes. The total transmission amplitude is given 
by the sum of the two amplitudes. 



t ^ ta+ tbC 



(A2) 



where $0 = h/e is the magnetic flux quantum. We as- 
sume that the influence of the qubit state on the ampli- 
tudes is the same for both junctions. Then, the ampli- 
tudes are related by a real constant c, ta.j = cttj for 
j — 0, 2. The total tunneling amplitude now becomes 
(up to an irrelevant global phase) 



t = v/l+c2(U + ta2(Tz) + taie-'^'x, (A3) 

where 5i = arg(l-|-ce~**/*''). Hence, this setup provides 
a way to obtain a tunneling Hamiltonian with tunable . 
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FIG. 3: (a) A finite pliase Si can be acliieved by using two 
junctions witli amplitudes ta{x,<Jz) and ti,{(Jz). The phase 
can be controlled by tuning the magnetic flux $. (b) A third 
junction tc and a two magnetic fluxes $1,2 yield control over 
both phases 5i,2- 



This has been important for the calculation of the noise. 
The setup can easily be extended to achieve the second 
tunable phase 62 which we used in the calculation of the 
current. Here, we need a third junction with an ampli- 
tude tc = tco, which is decoupled from both oscillator 
and qubit, and two magnetic fluxes $1,2- A schematic is 
shown in Fig. [5]d. 

Note that in order to observe Aharonov-Bohm phases, 
the phase coherence length (which can be on the order of 
microns) must exceed the size of the setup. Moreover, the 
change of the area due to the fluctuating position of the 
resonator is assumed to be negligible. Our analysis only 
covers the single- channel case. In the case of N chan- 
nels in the loop, the magnitude of the Aharonov-Bohm 
effect is reduced by a factor Finally, we assume 

that roundtrips of electrons are ruled out by the device 
geometry such that electrons leave the setup after pass- 
ing through either of the junctions. In this case, one can 
generally have \t{^)\^ ^ 



Appendix B: Current calculation 



The current operator is / 

A^tl^l^tpL), where A = jY and ^pa 



i^aiz = 0). The 
tunneling amplitude 7 is defined in Eq. (3). The av- 
erage current is in general time-dependent, and can be 
calculated by treating Ht as a perturbation to the non- 
interacting Hamiltonian Hq — Hei + Hr + Hq. In that 
case, the Kubo formula gives 



{i{t)) - -I / dt' {-^m [h^,{t) , i?T,i„t(i')] i*(o)), 

Jo 

(Bl) 

where the subscript "int" refers to the interaction picture 
and |vE'(0)) is the state of the total system at time t — 
0. We shall assume that |vE'(0)) is a product state of 
the fermionic ground state and the (possibly entangled) 
qubit-oscillator state. 

The current can also be expressed in terms of the state 
at time i, which to zeroth order in the tunneling ampli- 



in the tunneling amplitudes, the current becomes 

{I(t)) = 2 Re /" dr 
Jo 

If the time t is large compared to the timescales of the 
non- interacting system [1/^^, 57^^, (e^ + A^)^^/^], we can 
extend the r-integral to infinity. On the other hand, the 
above approximations are only valid for t smaller than 
any damping time scale, both due to the coupling to the 
leads (~ rj""^,r^^) and to other sources of dissipation 
not included in our model. 

Next, we define the electronic correlation functions 
Z^i-r,V) = {^lmL{~r)}{^PRmU-^)) 

ZU-T,V) = {^L{-r)i^lm{4'U~^)^Bm- (B3) 

The averages are taken with respect to the non- 
interacting ground state such that these correlation func- 
tions can easily be calculated. The chemical potentials of 
the left and right reservoirs differ by the applied voltage 
V . For simplicity, we focus on the simplest case of zero 
temperature and a constant density of states po- For zero 
voltage, one finds after Fourier transformation 

Z>{uj) = Z>{uj, V = Q)^ -2T:ple{-uj)uj 

Z< [uj) = Z<{uj,V = Q) = 27rpl0{uj)io. (B4) 

A finite voltage leads to additional phase factors 
Z>{t,V) = e-*^^Z>(r). By exploiting F^t = 1^ one 
then obtains the following expression for the average cur- 
rent 

(/(t)) =2Re |^°°dre^^^(z>(-r) (77^-^)) 

'Z<i-r){jH-r)l)), (B5) 

where the expectation value is in the state at time t. The 
time dependence of the operator 7 is in the interaction 
picture and is given by the free evolution of the oscillator 
and the qubit. 



P 

xcos{nt) H — ■ sin(r2t), 

mil 



e2 + A2 cos(2Ve2 + AH) 



Ae[l - cos(2Ve^ + AH)] 



A sin(2Ve2 + AH) 



A2 



(B6) 



The final expression (8) for the current can be ob- 
tained by using this, as well as [x,p] = 



tudes 7j is \^{t)) 



^(0)). Thus, to second order 



I, OiOj = 

6ij + iSijkCTk, and the assumption that the cut-off fre- 
quency of the density of states in the leads is much larger 
than V,n,Ve^ + A'^. 

In conclusion, given the state |^'(0)) at i = assumed 
above, the expression for the average current is valid for 

times v-\n-\ (e2 + A2)-i/2 « t < r^\r^\ 
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Appendix C: Noise calculation 

We now provide some details on the calculation of the 
symmetrised noise spectral density 



SM = l e-(({/(r),/(0)})), 



(CI) 



where ((/(r)/(0))) = (/(r)/(0)) - (/(r))(/(0)) and {•,•} 
denotes the anti-commutator. The time dependence is 
in the Heisenberg picture and the expectation values are 
taken in the state \^{t)) at time t. 



where as before /int(''') = e'^^"'^ I e '■^»'^. Applying the 
same approximations as in the calculation of the current, 
we find S{io,t) = {S{LJ,t) + S*{-uj,t))/2, where 



t)= dr e'(-+^)- (z> (-r) {^irh^) 

J — OC 

+ Z§i-r){lH-rh)). (C3) 



1. Perturbation theory 

To lowest order in the tunneling amplitudes, the noise 
spectral density can be approximated by 



S{u;,t) = - 



dr (f(t)|{/i„t(r),/(0)}|M/(t)), 

I 

(C2) 



S{uj, t) is an even function of cu. For \uj\, 17, A ^ \V\, it is 
frequency- independent and proportional to sgn(y)(/(i)), 
which is due to electron shot noise. Assuming e = and 
defining the function /(w, V) = + &{\V\ - \u\){\V\ - 
\ui\), the noise spectral density becomes 



S{uj,t) = npl 



l7o|(l7o| + l7i| cosSiix) + |72|cos52(a,))/(c<;,F) - 2n\ji\^xl + AA\j2\'' (a.) 
+ |7i I (|7o I cos 5i (x) + |7i I (a;') + |72 1 cos 621 (xa,)) (/(w, + 17) + /(w, V - fi)) 
+ 2a;^|7i| (|7o| smSi{p) - I72I sin52i(pa,)) (/(w, V + Q)- /(w, V - Q)) 
+ I72I (|7o| cos52(a.) + I72I + |7i| cos^2i(a;a^)) (/(w, V + 2A) + /{lo, V - 2A)) 
+ I72I (I70I sin (52 K) + |7i| sinS2i{xay)) (/(u;, ^ + 2A) - fico, V - 2A)) 



(C4) 



r 



We observe that S{uj,t) has kinks at |a;| = |y|, 
= |y±f2| and = |y±2A|, which can also provide 
information on the bipartite expectation value matrix of 
the oscillator-qubit system. 

As determining higher-order contributions to the noise 

spectral density is complicated using this approach, we 
now turn to an alternative method. 



sity matrix ps{t) = Trsp(t), will be determined by the 
Redfield equation. 



ps{t) = -i[Hs,ps{t)] 

/■OO 

- / dTTrB[HT,[HT{-T),ps{t)(3pB]]. (C5) 

Jo 



2. Born-Markov approximation 

In the case V » f2. A, e, the timescale of electron 
dynamics will be faster than the oscillator and qubit 
timescales. The electron reservoirs can then be regarded 
as a fermionic bath Hb = Hgi, coupled to the remaining 
system Hs = -\- Hq via the tunneling operator Ht- 
For small tunneling, the dynamics of the reduced den- 



The time-evolution of Ht is taken with respect to the 

uncoupled Hamiltonian Hg + Hg . Expectation values of 
system operators can in principle be calculated by using 
e.g. ^ {x{t)) = Trs[ps{t)x] and solving the resulting dif- 
ferential equation. Splitting the tunnel Hamiltonian into 
a system and a bath part, Ht = Aipl^ipii -|- h. c, we can 
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rewrite this as 

ps{t) = ~i[Hs.Ps{t)] 

^ dT[z>{-T,V)[A,A\-T)ps{t)] 

- Z<{-T, V)[A,ps{t)A\-T)] + h. c. }, (C6) 

where we used the bath correlation functions defined in 
Eq. (|B3P and the free evolution of the qubit-oscillator 
operators is given in Eq. (jB6[) . Therefore, A{t) contains 
terms oscillating at the frequencies and 2v? + A^. 
The r-integration in Eq. (IC6p can then be performed and 
its result can be absorbed into generalized transmission 
probabilities EJ^j where K and J run over the system 
operators, K,J £ {x,p,ax,cry,az}. The average current 
is defined as (/(t)) — T'i:s[ps{t)nR]- Using the Redfield 
equation, it can be written as 



W =-T.{^Kj{[nR,YK]Y^j) 

KJ 

-{YU[nR,YK]) +C.C.} 



KJ ' 



Y.{^KJ (KJ) - T^j (JK) +C.C.}, (C7) 



where we evaluated the expression using [Y, nn] = Y and 
yyt = yty = j -vvhile all other commutators of nji 
and Y with the system operators K, J vanish. The cal- 
culation of the transmission coefficients is a straightfor- 
ward task and the average current ior V > ft, y/e^ + A^ 
coincides exactly with the result from the perturbation 
theory. 



KJ 



In the following, we consider the case e — S2 = 0. 
The symmetrized frequency-dependent current noise is 
defined by Eq. (|Cip . In order to transform this into an 
expression which can be calculated using the Redfield 
equation, we use MacDonald's formula, see Eq. (11). The 
time-derivative of the cumulant of the number of trans- 
fered electrons ((n^)) = (nj^) — {njif' can be calculated 
using Eq. (IC6[) . The noise turns out to have a zero- 
frequency component proportional to the average current 
S{ljj = 0) = 2e (/) . After subtracting this term, one finds 



|((4,)) -(/> = ! (4) -2{nn) (/) -(/) 

/"p p T/p 

= W\J cos{5i)inRx)) + 2^UT^ns{n{5i)inRp)) + -^((n^.^^)) 

2 Ar 

+ -Y^{{nROx)) + ^^/vf^2{{nR(7,)) 

I v^r r / p p /p p ^ 

+ 2 J — ^—^ cos{5i)inBxaz)) - 2AW -i-^ s\w{5i){{nRX(jy)) + \ \/ sin((5i)((ni^pCT^)), (C8) 



where the oscillator coordinate was measured in units 
of its zero-point motion, x — x/xq^ p ~ 2xqp where 
xq = l/v2m?2. Next, we assume that the bare tunneling 
amplitude is large compared to the coupling to the qubit 
and the oscillator, ie. 70 ^ JiXo,j2- Since all of the 
above cumulants vanish in the absense of coupling to the 
qubit and the oscillator, they have to be on the order of 
these couplings 71,72- Therefore, cumulants containing 
quadratic operators in addition to ur are subleading and 
we can use 



dt 



((4)> - (/) = W^l^cosiSMnRx)) (C9) 



2yyfoI\nsm{6i){{nRp)) 
2AT2 



V 



)) + 4^vror2{{nR<jS. 



Since for a Gaussian distribution, all cumulants beyond 
quadratic order vanish, we call this the Gaussian ap- 
proximation. The next step is to calculate the time- 
dependence of the cumulants which can again be accom- 
plished using the Redfield equation (jC6|) . In order to sim- 
plify the expressions, we can assume (x) = (p) = (az) = 
since finite values for these expectation values would 
merely lead to a renormalized tunneling amplitude. As 
an example, consider the differential equation for {{urp)), 



j^iuRp)) = {■■■)- 2A^^ co^{5MnR<ry)) 



VT1T2 



sin((5i)((n_RCT^)) 



-^{{nRx))-T^{{nRp)), (CIO) 
where (• ■ ■ ) denotes a number of n/j-independent terms 
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and the initial condition reads {{njipfit=Q = 0. For the 
reason mentioned above, the terms containing cumulants 
of qubit operators are subleading and can be dropped. 
The same is true for the remaining differential equation, 
such that equations describing qubit cumulants and os- 
cillator cumulants decouple. This makes an analytical 



solution of the problem possible. The result is an expres- 
sion for the noise of the form S{uj) = 'S'x(w) {X) , 
where X denotes qubit or oscillator operators or prod- 
ucts of these. The prefactors Sx{'^) for the cross-terms 
turn out to be given by 



5,,^ = 2y2MrV\r2(A/l/) pa^iuj) + 2VP3{oj) cos((5i) + W^i^) cos(Ji) - WPi{uj) sin(5i)] 
S^,^ = 2^2MTxT2/V^ [%VToI5i{lo) cos(,5i) + Ta/^alw) sin(^i)] 

S,„^ = 2^J2MVT{C2 [-2r2(A/V^)2/33(w) cos((5i) + lWTof32{io) cos((5i) + 2Tona2{uj)\ 
Sp^^ = 2V2rori/Afr2(A/T/) [ai(c^) + 2/3i(cj) cos(5i) + I33{lu) sin((5i)] 
Sp^^ = ~^2Tir2/{VM)iA/V) [r2/33(w) cos(,5i) ~ 8yro/3i(w) sin((5i)] 
Spa^ = Ay/wf^r^To + 4/32 (w) sin((5i)] 

S^2,^ = 2Mrir2A/33(w) 

S,2^^ = Uly/ToT^TiA [-na2{io)sm{di) + 2V(3i{u)] 
S,2,^ = 4Afriv/VTV^[41//32(w) + ria2(t^)cos(Ji)] 

Sp2a^ = 

=2yr\^W(ri/Af)ai(t^)sin(5i) 
5s^^,„ - -2v/ror2/FriAai(c^)sin((5i) 

^Sx,a, = yfiirWri [2Vai{ij) cos(Ji) + l]a2(t^) sin(<5i)] , (Cll) 



where the functions Q!i,2(w) and /3i,2,3('^) were defined in Eq. (12). 
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